Mon. Not. R. Astron. Soc. OOP. [TlfTTI C2009') Printed 12 April 2009 (MN I^TeK style file v2.2) 



Determination of the intrinsic velocity field in the M87 jet 
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ABSTRACT 

A new method to estimate the Doppler beaming factor of relativistic large-scale jet 
regions is presented. It is based on multiwaveband fitting to radio-to-X-ray continua 
with synchrotron spectrum models. Combining our method with available observa- 
tional data of proper motions, we derive the intrinsic velocity as well as the viewing 
angles to the line of sight for eight knotty regions down the M87 jet. The results fa- 
vor the 'modest beaming' scenario along the jet, with Doppler factors varying between 
~ 2 — 5. The inner jet of M87 suffers sharp deceleration, and the intrinsic speed remains 
roughly constant down the outer jet. The orientation of the inner jet regions is fully 
consistent with the result of 10°-19°to the line of sight suggested by previous Hubble 
Space Telescope (HST) proper motion studies of the M87 jet. The outer jet, however, 
shows systematic deflection off the inner jet to much smaller inclination {9 ^ 10°). 
Further calculation of knot A suggests this deflection can be regarded as evidence that 
the outer jet suffers some departure from equipartition. The nucleus region of the M87 
jet should have a viewing angle close to its first knot HST-1, i.e. 9 ~ 15°, which favors 
the idea that M87 may be a misaligned blazar. This work provides some hints about 
the overall dynamics of this famous extragalactic jet. 

Key words: galaxies: active - galaxies: individual (M87) ~ galaxies: kinematics and 
dynamics - galaxies: jets. 



1 INTRODUCTION 

Since the discovery of the first extragalacti c large-scale jet 
in a nearby active galaxy M 87 (|Curtislll918l ). intense efforts 
are made to explore the underlying physics of the relativis- 
tic jets during the past decades. The related topics concern 
about the matter content, particle acceleration, coUimation 
mechanism and magnetic field configuration inside the jet 
flow. Among these problems, one of the fundamental issues 
is the intrinsic velocity field of the jet. Lack of an effective 
method to determine the speed and orientation of the suc- 
cessive regions in large-scale jet, it is difiicult to constrain 
further the mechanism on overall dynamics of the bulk jet 
flow. 

In this work, a new method is presented to estimate the 
Doppler beaming factor in a region of large-scale jet. This 
method is based on the spectral flt of the multi-waveband 
continua with standard synchrotron radiation models. Com- 
bining the data of proper motions, we readily derive the 
intrinsic bulk velocity as well as the angle to line-of-sight 
of the jet flow. Making application of this new method, we 
successfully obtain the intrinsic velocity field in M 87 jet. 
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Our result agrees well with that of the n ew variability mon i- 
toring program toward the same object (jHarris et al.ll20()6l l. 
which indicates the reliability of our approach. The method 
is simple and only requires non-simultaneous spectral energy 
distribution (SED) data sets, therefore is hopeful to be of 
general interest for those synchrotron-dominated jets. 

For the rest of this paper, we firstly make synchrotron 
spectrum model fitting to radio-to-X-ray SED of M87 jet 
in §2, then introduce the new method and derive the dis- 
tribution of Doppler beaming factor along M 87 jet in §3. 
Coupling the data of proper motions, we obtain the intrin- 
sic velocity field of M87 jet in §4. Related discussions are 
also presented in this section. Finally, the main results are 
summarized in 55. 



2 SYNCHROTRON SPECTRUM MODEL 
FITTING TO M87 JET 

2.1 The Synchrotron Radiation Models 

The high spatial resolution of Chandra X-ray observatory 
(FWHM = 0.5") has opened a new era on the research of 
extragalactic large-scale jets at high energy waveband. At 
the time of preparing this paper, 91 radio-loud AGNs are re- 
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ported detection of X-ray counterparts of radio jets on large- 
scales (see also |http : //hea-www . harvard . edu/X JET/ ) . The 
number of detected X-ray jets increases rapidly, suggesting 
that X-ray emission from jets is a common feature in radio 
galajdes and quasars. If combined with observational data 
in radio and optical band (sometimes also with available in- 
frared and ultraviolet data), we can construct broadband 
SED and carry out multi- waveband modelling of the jet re- 
gions. Such a procedure may provide us with useful infor- 
mation of jet physics. 

Polarization observations already confirm the syn- 
chrotron nature of t he radio-optical emi ssions from ex- 
tragalactic jets (e.g., iPerlman et al.l Il999l for M 87 jet). 
The X-ray emission process of large-scale jets, however, re- 
mains as an open question. For nearby FR I radio galax- 
ies, the radio-to-X-ray SED from the large-scale knots is 
consistent with a single smoothly broken power-law spec- 
trum, and the expected inverse-Compton scattering (IC) 
fluxes of seed photons (jet synchrotron photons or cosmic 
microwave background photons) always underpredict the 
observed X-ray fluxes. This implies the synchrotron ori- 
gin of the X-ray em issio ns observed in FR I jets (e.g.. 



MarshaU et al.1 12002| and IWflson fc Yand I2OO2I for M 87, 
Hardcastle. Birkinshaw fc Worralll I2OOII for 3C 66B). The 



recent statistical research of a large X-ray jet sample 
iKataoka &: Stawar j|2005l ') also favors this physical picture. 
On the other hand, the X-ray emissions from large-scale 
quasar and FR II jets, however, are usually interpreted as 
inverse-Compton scattering of CMB photons within a highly 
relativistic jet ( Harris fc Krawczynski 2002). 

As discussed bv ISchwart3 (|2002l '). the IC/CMB model 
for the X-ray emissions of large-scale quasar jets predicts 
two main differences between nearby and distant objects. 
Due to the enhancement of the CMB energy density with 
redshift by a factor of (1 -I- 2;)^, the ratio of X-ray to radio 
fluxes of the resolved parts of the quasar jets should increase 
in distant sources compared to nearby ones. The other effect 
is that we should observe some systematic flattening of the 
optical-to- X-ray spectral index (aox) in distant unresolved 
quasar cores compared to low-redshift quasars, because of 
the increasing contribution from the unresolved portions of 
the jet to the X-ray flux of the quasar core. 

We not ice, however, the recent observational efforts of 
iLopez et al.l |200Q) to search the latter effect. They carry out 
a Chandra snapshot survey of representative high-redshift 
radio-loud quasars [z ~ 3.5 — 4.7) selected from the Parkes- 
MIT-NRAO sample. The survey does not detect any system- 
atic flattening of the optical-to-X-ray spectral index (qqx) 
of the unresolved quasa r cores compared t o low-redshift 
quasars. The results of iLopez et al.l (|2006l ) suggest that 
kiloparsec-scale X-ray emission of quasar jet is not dom- 
inated by inverse Compton scattering of CMB seed pho- 
tons off jet electrons. Meanwhile, the IC/CMB model also 
meets difficulty in exp laining the X-ray emissions from sev- 
eral FR II jets (e.g., iKraft et al.1 12005| for 3C 403). It is 
also worthy to notice the most recent discovery of the X-ray 
counterjet in the FR II source 3C 353 detected by Chandra 
iKataoka et al.ll2008l ) . The research of lKataoka et al.1 (|2008l ) 
suggests that this detection is inconsistent with the IC/CMB 
model, and instead implies a synchrotron origin of the X-ray 
jet photons. Thus, even for powerful quasar/FR II jets, the 



synchrotron radiation mechanism for X-ray photons keeps 
as a plausible interpretation. 

In general, synchrotron origin for X-ray emission of 
large-scale jet faces much more strict challenge in particle 
acceleration than IC/CMB model does. Firstly, the particle 
acceleration in synchrotron X-ray jets has to be fast enough 
to generate ultra-relativistic electrons with maximum ener- 
gies ~ 10 — 100 TeV (i.e., one to two orders of magnitude 
higher than that required in the IC/CMB model). Secondly, 
the electron acceleration processes in synchrotron X-ray jets 
have to operate continuously within the whole volumes of 
the jets. Such a strict requirement that synchrotron X-ray 
jets should work as powerful particle accelerators is naturally 
supported by the idea that large-scale jets are very hopeful 
sites for yielding Ultra High Energy Cosmic Ray particles 
(UHECRs, see e.g. [c asse fc Marcowit"hll2005l and references 
therein). 

Based on the observations mentioned above, we choose 
synchrotron radiation model for our multi-waveband re- 
search. Currently, there are three standard synchrotron 
spectrum models available for statistical flt of SED: (1) 
The Jaffc & Pcrola (1973) model (hereafter JP), which as- 
sumes the pitch angle distribution of relativistic electrons is 
continuously isotropized after an initial injection of single 
power-law electron energy distributio n. (2 ) The Kardashev- 
Pacholczyk model (jKardashevI Il962l and iPacholczvkl Il970l . 
hereafter KP), which allows only the evolution of high en- 
ergy tail of electron energy distribution following an ini- 
tial single-power law electron injection. Since no pitch-angle 
scattering of the radiating electrons is allowed, it is unlikely 
to be the real case. This model, however, can give good 
fits to the observational data in many cases, therefore is 
relatively common in the je t research community. (3) The 
Continuo us Injection model (|Heavens fc Meisenheimeij[l987l 
and Mcis enheimer et al.l Il989l . hereafter CP), which allows 
continuous electron injection with single power-law energy 
distribution into the emission region. 

The spatial extent of the X-ray emitting regions in 
some FR I jets resolved by Chandra implies that continu- 
ous acceleration and injection of high-energy electrons may 
oc cur in these objects. On the other hand, recent work 
of iPerlman fc WilsonI (|2005h on M87 jet suggests that JP 
model underpredicts the X-ray fiux by many orders of mag- 
nitude and the theoretical slope at X-ray energies is much 
larger than those observed. Therefore we do not consider the 
JP model in our work. All of the multi-band fits are made 
using KP and CI models. For all of the synchrotron spec- 
trum models above, the dominated cooling process of high- 
energy electrons responsible for X-ray emissions should be 
synchrotron radiative energy loss. 



2.2 Result of Multi-waveband Fit 

M 87 is one of the nearby FR I radio galaxies (distance = 
16 Mpc, with a scale of 1"= 77.6 pc). The rich observational 
data of its famous jet enable us to do further multi-band fit. 
The program of synchrotron spectrum fit wo use is written 
by C. Carilh and J. P. Leahy (Carilh et al. 1991.; , Leahy 199ll). 
The flt is done in two steps. In step 1, we get an initial guess 
of the position of th e break frequency by emp loying a Myers 
and Spangler test (|Mvers fc Spangleil 19851 ) on the multi- 
waveband data sets. Then in step 2, a Marquardt non-linear 
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test is done numerically (jPress et al.lll987l ) to reach the 
best fitting of the SED data. The tests are model dependent. 
The electrons responsible for higher frequency synchrotron 
radiation suffer more rapid synchrotron loss, thus there is 
a spectra steepening toward the high frequency band. For 
CI model, there is a break in spectral index of Aa = 0.5 
between the low and high frequency spectra. In the case 
of KP model, Aa depends on the index of electron energy 
distribution. 

In Figure[T] we plot results of CI as well as KP model fit- 
ting for the radio-to-X-ray SED of nucleus and eight knotty 
regions. It is worth to notice that knot HST-1 has started 
to flare (in radio, o ptical and X-ray waveband) around 2002 
l|Harris et al.ll200^ . Since we aim to do 'quiescent state fit' 
of the non-simultaneous SED, the broad-band data for this 
knotty region of M87 jet are carefully selected to be prior 
to the flare. The VL A observational data point at 15 GHz 
is from [Zhoul (|l998h . and optical -near-infrared data p oints 
are obtained with HST on 1998 jPerlman et al.ll20qial). A s 
to the X-ray spectrum, it is from iPerlman fc WilsonI l|2005[ ) , 
which is the reanalysis o f early Chandra observations of M87 
jet taken on 2000 July (| Wilson fc Y ang 2002). The observa- 
tions above loc ate prior to the outbu rst of HST-1 knot. 

Moreover, IWaters fc Zepj (|2005l ) present new ultravio- 
let HST observations of M87 jet taken on 2001 February 23. 
Their results suggest that knot HST-1 reveals a significant 
increase in the brigh tness from late 2001 to the present. 
IWaters fc Zepj (120051 ) try to do best fit to the observed 
radio-to-X-ray SED of HST-1 with standard synchrotron 
spectrum models, and found the new UV point exceeds the 
fit by a large amount independent of the models. Their work 
implies HST-1 knot has started its early evolution of the flare 
by the time of their observations. Therefore, we include the 
UV observational data points of IWaters fc Zepl l|2005l ) for 
most M87 jet knots (no significant variability found up to 
date) except HST-1. 

Table [1] lists the best-fit parameters for synchrotron 
spectrum models. During the fitting procedure, we set the 
input low frequency spectral index Oinj as the observed radio 
through optical spect ral index a^o (the same treatment as 
iPerlman et al.ll2001al ). Here the convention of synchrotron 
spectrum is Si, oc v'". 

In Figure [S] we plot the derived break frequency vb 
along M87 jet. Our result reveals a general trend of decline 
on the break frequency of la r ge-sca le knots, similar to that 
discovered bv lMarshall et al.l 1 2002f) . We find, however, a sig- 
nificant increase of vb from jet nucleus to first knot, HST-1, 
by both of the model fit t ing. W ith CI model fitting to mid- 
IR SED. IPerlman et"al] l|2001bf ) derived break frequency of 
M87 nucleus, vb = 2.8 x lO^^Hz. The good agreement be- 
tween our derived ub (2.22 x 10^^ Hz) of the nucleus based on 
multi-waveband fitting and that of IPerlman et"al] l|2001bl ') 
suggests that the increase-decline trend is reliable. 

Here we discuss further the origin of broad-band spectra 
of the unresolve d nucleus of M8 7 jet. T here is some debate 
on this problem. [Marshall et ahl l|2002l ) present high resolu- 
tion X-ray image of M87 jet using Chandra. They find that 
the core fiux is significantly larger than expected from an 
advection accretion fiow and the spectrum is much steeper. 
Their results indicate that the X-ray emission of the M87 
nuc leus is due to synchro tron radiation from a small scale 
iet. lWilson fc Yang l|2002l ) also give Chandra X-ray imaging 



and spectroscopy of M87 jet. In view of the similar spectra 
of the nucleus and jet knots, and the high X - ray flu x of the 
knots closest to the nucleus, IWilson fc Yand (|2002l ) suggest 
that the X-ray emission of M87 nucleus may actually origi- 
nate from the pc- or sub-pc-s cale jet rather than the accre- 
tion disc. On the other hand, iDi Matteo et al.] (|2003t ) argue 
that the observed radio-to-X-ray spectrum of M87 nucleus 
is consistent with that predicted by an radiatively-inefficient 
accretion disc. However, they can not also rule out the pos- 
sibility that the X-ray emission of the nucleus is dominated 
by jet emission. 

Considering this problem, it i s also worthy to notice 
the latest resear ch of iLenain et al. | l|2008l '). Similar to pre- 
vious researches, iLenain et al.l ( 20081 ) suggest that standard 
one-zone jet model can describe roughly the radio to X- 
ray SED of M87 nucleus (see also figure(la) of ILenain et al.l 
(2008)). They find, however, the single-zone approach can- 
not describe correctly the very high energy (VHE) emis- 
sion of M87 detected by High Energy Stereoscopic System 
(HESS) telescope array. Their further research indicates that 
a multi-blob jet model of M87 nucleus may do better at TeV 
energy band. We note, however, there is some controversy on 
the origin of the unresolved TeV 7- ray source found in M87 . 
For example, detailed modelling of iHonda fc Hon"dal (|2007| ) 
suggests that the VHE emissions may originate from bright 
knot A of M87 jet. Therefore, we argue that a single-zone 
approach of jet model is tolerable if our research interests 
focus only on radio to X-ray SED of M87 nucleus. 

If the broad-band emissions of M87 nucleus mainly 
come from the unresolved compact jet region, it is nec- 
essary also to evaluate the influence of synchrotron-self- 
absorption effects at radio frequencies upon the result of 
synchrotron model flts. We flnd out that the radio spectrum 
of M87 nucleus is roughly consistent with a single power- 
law, i.e. Si, oc v'~'^'^ . The flat radio spectrum implies that 
M87 nucleus is only partially opaque at the observed ra- 
dio frequencies, and is far beyond the optically t hick regime 
of sy nchrotron-self-absorption, i.e. Sv oc v'^'^ (jPacholczvkl 
I1970I ). Consequently, the synchrotron model fit to radio- 
to-X-ray SED neglecting the synchrotron-self- absorption ef- 
fect should alter the derived break frequency only slightly, 
e.g. our derived i^B,ci = 2.22 X lO^^Hz is i n good consis- 
t ent wi th fB.ci = 2.8 X lO^^Hz derived by jPerlman et al.l 
(2001b) based only upon mid-IR SED of M87 nucleus. We 
will further indicate that the deduced Dopper beaming fac- 
tor is insensitive to the exact value of break frequency (see 
§3), therefore the simplification of broad-band fit to M87 
nucleus at the radio frequencies should not alter the result 
significantly. 

Inspection of Figure [1] and Figure [2] also indicates that 
there is systematic difference between CI and KP model fit- 
ting. KP model predicts higher break frequency than CI 
model does. Throughout M87 jet, CI model overpredicts 
the X-ray ffux and the theoretical X-ray spectrum is flat- 
ter than that observed by Chandra. Vice versa, KP model 
underpredicts the X-ray flux, and the predicted X-ray spec- 
trum is steep er than tha t observed. The same result is also 
achieved bv iPerlman fc Wilson (.2005il . The systematic de- 
viation of the two models off the X-ray spectra implies that 
the real position of vb may deviate also from standard CI 
and KP cases, or a modified non-standard synchrotron spec- 
trum model should work better. 
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Figu re 1. CI (so lid li ne) and KP (dashed line) synchrotron s pectrum model fits for nine regions alon g M87 jet. Data of knots are from: 
VLA jZhoulll998l: HSTjPerlman et"alll2001al . [Waters fc Zepj |2005: Cha ndra-IPerlman fc WilsonI l2005l . Data of nucleus are from: VLA, 
HST- ISparks. Biretta fc Macchettolll996l and references therein; Ciiandra JPerlman fc Wilsonll2005[ 



Table 1. Best-Fit Parameters for Synchrotron Spectrum Models 



Jet Region 




CI Model 








KP Model 






Nucleus 


0.70 




= 2.22 ±0 


12 


X 


IQi^Hz 


3.29 ± 0.07 


X 


IQISHz 








= 119.47 








56.05 






HST-1 


0.71 




= 1.84 ±0 


09 


X 


IQi^Hz 


5.91 ± 0.19 


X 


10"Hz 








= 7.51 








5.69 






D-East 


0.70 


I'B 


= 2.58 ±0 


10 


X 


IQi^Hz 


6.44 ± 0.08 


X 


IQi^Hz 






X^ 


= 27.91 








36.30 






E 


0.71 


I'B 


= 4.36 ±0 


22 


X 


10i'""'Hz 


8.17 ±0.14 


X 


lO^'^Hz 






X^ 


= 15.26 








20.83 






F 


0.69 


I^B 


= 8.95 ±0 


43 


X 


IQi^Hz 


1.95 ± 0.04 


X 


lO^'^Hz 






X^ 


= 197.93 








14.41 






A 


0.67 


I^B 


= 2.16 ±0 


06 


X 


IQi^Hz 


1.42 ± 0.02 


X 








X^ 


= 610.35 








18.81 






B 


0.67 


I'B 


= 2.38 ±0 


07 


X 


10"Hz 


5.45 ± 0.10 


X 


10i"^Hz 






X^ 


= 764.86 








8.58 






C-1 


0.69 


I^B 


= 2.42 ±0 


06 


X 


IQl^Hz 


3.30 ±0.04 


X 


lO^^Hz 






X^ 


= 2554.64 








59.76 






C-2 


0.68 


I^B 


= 9.64 ±0 


28 


X 


lOi^'Hz 


1.34 ±0.01 


X 


lO^f^Hz 






X^ 


= 508.73 








92.74 







In order to explain such disagreement between the pre- 
dicted and observed X-ray spectrum in M87 jet knots, sev- 
eral possible solutions to this problem are proposed. For 
example, if considering synchrotron cooling as well as elec- 
tron acceleration processes in a non-uniform magnetic field, 
iBickneU fc BegelmanI {l99&) suggest it is possible to achieve 
a larger break in spectral index than sta ndard CI model, 
i.e. Aa > 0.5. With a similar scenario, iHonda fc Hondal 



obtain a modified X-ray spectral index consistent 
with Chandra observations of M87 knot A. 

During th e preparation of th is paper, we also notice the 
latest result of lLiu fc ShenI (j2007l ). They propose a modified 
CI model to explain the radio-to-X-ray continua in six knots 
of M87 jet. Considering the thin acceleration region (i.e. 
shock front) locating at the immediately upstream of the 
main emission blob, the broadband spectra can be fit much 
better than standard CI model. 
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Figure 2. Derived break frequency through CI and KP model fitting versus the projected distance from jet nucleus. 



3 DERIVING DOPPLER BEAMING FACTOR 
ALONG M87 JET 

Through synchrotron spectrum model fitting, we can esti- 
mate brealc frequency (vb) of the observed SED. This fre- 
quency provides us with important information of the jet re- 
gion. In general, we consider injection of a power-law distri- 
bution of relativistic electrons into an emission region dom- 
inated by synchrotron cooling process. Then the observed 
break frequency will correspond to a characteristic break 
energy of electron population (7 = 7br), where synchrotron 
radiation loss is ba lanced by escape of electrons from the 
emiss ion blob (e.g. Ilnoue fc Takaharal Il996l . iKataoka et al.l 
I2OOOI ). 

Therefore we have an relation at the break frequency 



isyn(l'B) = tc 



(1) 



where tsyn(i^B) is synchrotron cooling time at the break fre- 
quency, tcsc is diffusive escape time scale of electrons from 
the emission region. 

The synchrotron cooling time can be further expressed 



tsyn(VB) = 5.08 X lO^^B 



(2) 



where the magnetic field B, is in ^G, and the break fre- 
quency i/B, is in GHz. 

Up to date, we know little about the random motion of 
electrons moving in the emission blob. Thus the escape time 
of electrons, tcsc, is quite uncertain. In general, it should be 
longer than hght travel time over the source {R/c, where 
R is the size of the emission region and c is speed of light). 
The escape time, however, is unlikely to be much longer than 
R/c. Otherwise, no spectral break (ub) will be observed in 



the SED, which is not the fact o f the o bservations. We then 
adopt the treatment of iKataokal (| 19991 ) as 



tcsc = vR/c, (1 < f? < 10) 



(3) 



where is a dimensionless parameter. To estimate the influ- 
ence of escape time upon the uncertainties of Doppler factor, 
we further assume that tesc is Gaussian distributed in this 
regime (mean value 5.5R/c, ±3(t range), i.e. 77 = 5.5 ± 1.5. 

The equations above are all written in jet frame. 
Considering the relativistic bulk motion of the jet re- 
gion described with a certain Doppler beaming factor, 5, 
there is a Lorentz transformation of the physical quan- 
tities from the source (jet) frame to the observer frame. 
According to the formulae p resented in appendix es of 
Harris fc Krawczvnskil ||200^ ) , IStawarz et all (120031 ) and 
Begelman. Blandford fc Rees (|l984 ). we substitute ^ and 
((31 into (HI and rewrite the physical quantities in the ob- 
server frame. Thus Doppler beaming factor appears on both 
sides of the relation, which naturally leads to the final for- 
mula of S (see Appendix A for details) 



vRob 



1.61 X 1021 



7/18 



(4) 



where -Robs is the mean observed radius of emission region 
in cm, c is speed of light in cm • s~ , I'B.obs is the observed 
break frequency in Hz, B{\) is the equipartition magnetic 
field in /iG calculated for no beami ng {5 = 1). We com - 
pute -8(1) using eq. (A6) of .Harris fc Krawczvnskil (|2002l ). 
which requires only total synchrotron luminosity and -Robs of 
the source. We can fit the radio-to-X-ray spectrum approx- 
imately with a single power law to estimate the observed 
bolometric luminosity of the jet region, following the treat- 
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ment of [Harris etU] (120031 ). On the other hand, if we have 
FWHM measurement of the major axis (&a) and minor axis 
{Oh) for a jet region from high resolution imaging obser- 
vations, we can use the quantity \/9a, x 9^, to estimate the 
mean angular diameter of the region, which in turn gives 
^obs with the known distance of the object. Due to the ex- 
treme proximity of M87 radio galaxy (distance = 16 Mpc), 
we neglect the redshift correction in equation Q, which re- 
quires an additional factor {1 + z) multiplied at the right 
side of that equation. 

Although KP model gives better fits to broad-band 
spectra of M87 jet, this scenario is somewhat unrealistic. 
The spatial extent of the X-ray emitting regions in FR I 
jets resolved by Chandra implies that continuous acceler- 
ation/injection of high-energy electrons may occur in such 
objects. Therefore, the break frequency used in equation (|4]) 
to estimate Doppler factor sho uld be derived from a more re- 
alistic modified CI r nodel (e.g..lLiu fc Shenl2007l ). According 
to the definition of iLiu fc Shen 1 20071 ). their second break 
frequency, pb2, corresponds to the physical break frequency 
which relates to the break energy in the spectral energy dis- 
tribution of relativistic electrons. They derive 1^32 for five 
knotty regions (D, F, A, B, CI), which agree with our i/b,kp 
around a factor of 2. Furthermore, they also report the same 
decline of that break frequency down M87 jet. 

On the other hand, we find in equation Q, 5 ~ z/g 
which suggests the value of Doppl er factor is insensit ive to 
exact position of break frequency. iLiu fc ShenI l|2007l ') only 
give physical break frequency for five knots of M87 jet. To 
keep a standard treatment, we use i^b.kp to approximate the 
physical break frequency derived from modified CI model in 
all nine regions along M87 jet. 

We list in Table [5] the related parameters as well as de- 
rived Doppler factors for nine successive regions down M87 
jet. Here i?obs is estimated us ing Chandra data of FWHM 
component sizes (see Table 1 of lPerlman fc Wilson 2005), for 
a scale of 1"=77.6 pc. We notice that iKataoka fc Stawarz, 



//5 



(|2005l ) calculate the equipartition magnetic field for three 
regions of M87 jet, i.e., knot HST-1, A and D. Their de- 
rived values of B{1) which is based only upon chromatic 
radio luminosity, agree with our result within a factor of 2. 
Therefore our derived parameters are reliable. 

The related parameters are substituted into Q, then 
the values of Doppler beaming factors are derived and listed 
in last column of Tabled The uncertainties of 5 are calcu- 
lated using standard theory of error estimation. Our result 
provi des strong support to the 'modest beaming' scenario 
(e.g., [Harris et al.l 12003 ) in this typica l FR I radio galaxy . 
It is worthy to note the new result of [Harris et al.l (20061 ). 
They give the primary observational result of the long-term 
Chandra monitori ng program of the variable knot HST-1. 
[Harris et al.l (|2006l ) report the X-ray decay time of HST-1 for 
the intensity to drop by a factor of 2 range from 0.2 to 0.3 yr. 
Their detailed analysis suggests that small values of Doppler 
fa ctor, i.e., 5 ~ 3 — 5 , are favored for this knot. The result 
of [Harris et al.l (|2006[ ) is well consistent with our indepen- 
dent resuh for knot HST-1, i.e., S 3.0 6 - 4.08. Moreover, 
[Podson. Edwards fc HirabavashU (|2006[ ) carry out detailed 
VSOP observations toward the nucleus region of M87 jet. 
Their study suggests that the M87 nucleus is not strongly 
Doppler boosted, with a Doppler beaming factor somewhat 
larger than 1.6. Our derived S ~ 1.79—2.95 for nucleus region 



agree with their work again for this region. The consistency 
above then provides strong support to our new method. 



4 DETERMINATION OF THE INTRINSIC 
VELOCITY FIELD AND DISCUSSIONS 

Once S and the observed apparent speed /3app of a jet region 
are available, we can compute the bulk Lorentz factor F and 
angle to line-of-sight 9 using eq. (10) and (11) of [Piner et al.[ 
(2003) as following: 



r = 



25 



and 



arctan 



2/3a 



/3|pp + 



(5) 



(6) 



M87 owns detailed observations of proper motion for 
most of its jet knots. In Table O we list the observed data 
of apparent velocity /Sapp. Coupling the already derived 
Doppler factor 5, the bulk Lorentz factor and viewing angle 
of the knotty regions can be calculated with (O and © (see 
also Table [3} . The overall results are visualized in Figure O 

First inspection of Table [3] indicates that the viewing 
angles of the inner jet regions are fully consistent with 10°- 
19°to line-of-sight sug gested by previous HST proper mo - 
tions study of M87 jet l|Biretta. Sparks fc Macchetto[[l999l ). 
For HST-1 knot, we use proper motion data up to ~ 6c, 
which are ba sed on yearly HST optical monitoring from 
1994 to 1998 ([Biretta. Sparks fc M acchctto 1999). Recently, 
[Cheung et al. 1 (|2007[ ) carry out VLBA monitoring of the 
same jet region from 2005 to 2006, and report apparent ve- 
locity up to (4.3 ±0.7)c, somewhat lower than previo us HST 
measurement of [Biretta. Sparks fc Macchettd ([l999l ). Com- 
bine this new apparent velocity with our derived Doppler 
factor of HST-1 knot, we find V = 4.51 ± 0.85 and 6 = 
15.9 ± 1. 9 deg. Our result ag rees well with F 4.4 sug- 
gested bv [Cheung et al.[ (|2007[ ). A comparison between two 
sets of results of HST-1 knot in two monitoring intervals sug- 
gests that the orientation of this jet region remains roughly 
the same, and the intrinsic velocity of jet flow suffers some 
decrease. Such a decrease of bulk velocity of the first knot 
of M87 jet may imply some reduced activity of its central 
engine. 

As alrea d y not ed by many researchers (e.g., 
[Perlman et al.l [2001al ). the morphology of the outer 
M87 jet (knot A, B, C-1, C-2 in Table (Sjl is quite different 
from that of the inner jet (knot HST-1, D-East, E, F). 
Inspection of Table [3] and Figure [3] indicates that there 
is systematic difference on the intrinsic velocity as well 
as orientation between these two segment of jet flow. In 
Fig. ((Sja), we find the overall distribution of Doppler factor 
down M87 jet is quite flat, supporting the 'modest beaming' 
scenario. The inner jet, however, suffers sharp deceleration. 
On the other side, the intrinsic velocity of the bulk jet flow 
remains roughly constant in the outer jet (see Fig. (|3];)). 
In fact, the similar general trend of /3app distribution (see 
Fig. ([SK)) provides us with some hints. 

In general, the early deceleration occurred in typical 
FR I jets such as M87 jet is believed to be caused by contin- 
uous mass loading through entrainm ent of arnbient gaseous 
medium during jet propagation (e.g. [BicknellllT993 ). Based 
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Table 2. Related Quantities and Derived Doppler Factors Down M87 Jet 



Jet Region 


(arcsec) 


6b (arcsec) 


^f-obs(pc) 


S(1)(mG) 


5 


Nucleus 


0.70 ±0.04 


0.58 ±0.03 


24.7 ± 1.0 


699 ± 264 


2.37 ±0.58 


HST-1 


0.79 ± 0.06 


0.60 ±0.04 


26.7 ± 1.4 


237 ± 38 


3.57 ±0.51 


D-East 


0.85 ± 0.06 


0.73 ±0.05 


30.6 ± 1.5 


285 ± 121 


2.72 ± 0.74 


E 


1.4 ±0.2 


0.9 ±0.1 


43.6 ±3.9 


183 ± 81 


2.53 ±0.71 


F 


1.6 ±0.2 


1.2 ±0.2 


53.8 ± 5.6 


211 ± 143 


2.26 ±0.93 


A 


0.98 ± 0.04 


0.91 ±0.04 


36.6 ± 1.1 


547 ± 362 


3.18 ± 1.27 


B 


1.7±0.2 


1.3 ±0.2 


57.7 ±5.6 


319 ± 250 


2.30 ± 1.08 


C-1 


0.49 ± 0.00 


0.49 ± 0.00 


19.0 ±0.0 


871 ± 668 


2.44 ± 1.12 


C-2 


0.49 ± 0.00 


0.49 ±0.00 


19.0 ±0.0 


556 ± 355 


2.46 ± 0.95 



9a and 6b a^re major an d minor component sizes (FWHM), respectively, from 
Chandra observations l lPerlman fc Wilsonll2005l) . Robs is mean observed ra- 
dius of emission region, B{1) is equipartition magnetic field calculated for no 
beaming, i.e. Doppler beaming factor 5 = 1. 



Table 3. Derived Intrinsic Velocity and Orientation of M87 Jet 



Jet Region 


/3app 


Reference 


r 


e (deg) 


HST-1 


6.14 ±0.58 


B99 


7.21 ± 1.12 


13.9 ± 1.2 


D-East 


2.99 ±0.28 


B99 


3.19 ±0.33 


21.3 ±5.1 


E 


3.92 ±0.80 


B99 


4.50 ±1.36 


20.7 ±3.8 


F 


0.86 ±0.23 


B95 


1.51 ±0.32 


19.5 ± 16.0 


A 


0.5090 ±0.0015 


B95 


1.79 ±0.56 


6.2 ±5.3 


B 


0.62 ±0.05 


B95 


1.45 ±0.40 


14.9 ± 15.1 


C-1 


0.11 ±0.04 


B95 


1.43 ±0.46 


2.54 ± 2.94 


C-2 


0.11 ±0.04 


B95 


1.44 ±0.40 


2.49 ± 2.46 



/3app is observed apparent speed in unit c, P and 6 refer to bulk Lorentz 
factor and viewing angle to line-of-sight of jet flow, respectively. Reference s 
of pr oper motion observations are : B99 llBiretta. Sparks fc Macchetto|[l999l) . 
B95 llBiretta. Zhou fc Owenlll995h . 



upon the derived velocity field here, we will explore the dy- 
namics of deceleration of M87 jet deeply in a later paper (C. 
Wang et al., in preparation) that employs full treatment of 
special relativity hydrodynamics. 

The orientation of inner jet and outer jet reveals differ- 
ent pattern also (see last column of Table [3] and Fig. ((HJi) . 
For the inner jet (knot HST-1, D-East, E, F), our de- 
rived viewing angles are fully consistent with 10°-19°to 
line-of-sig ht suggested by the HST proper motion study of 
iBiretta. Sp arks fc Macchet tol 1(199^ in their Table 3. As to 
the outer jet (knot A, B, C-1, C-2), however, we find sys- 
tematic defiection to much smaller viewing angles( Despite 
of the large uncertainties of viewing angle of knot B, we in- 
fer its most plausible orientation lies intermediately between 
that of knot A and knot C-1 considering smooth propagation 
of jet fiow). 

Lack of exact m easurement of proper motion (e.g., 
ILv. Walker, fc Junoij |2007) give VLBI detected apparent 
speed of core jet as (0.25 — 0.40)c, while they argue it should 
be treated only as lower limits), we do not calculate the 
intrinsic velocity and orientation for the nucleus region of 
M87 jet. We argue, however, the viewing angle of M87 nu- 
cleus should be very close to the first knot, HST-1, i.e. 
d ~ 15° based upon mean estimation of values of 61 ~ 14° and 
16° desc ribed above. Using t he observed TeV band spectra 
of M87, iLenain et al.l (|2008l ) also constrain the orientation 
of its core jet as 6^ ~ 15° . Their results present good support 



to our viewpoint. Therefore, the idea that M87 may be a 
misaligned blazar is favored by our research. 

Here we explore further the somewhat strange orienta- 
tion of the outer jet. We obtain small viewing angles (^ 10°) 
for the outer parts of M87 jet. The results confiict, however, 
with much larger inclination sugge sted by previous research 
(e.g., iBicknell fc BegelmanI l| 19961 ) indicate the orientation 
of knots A and B should range between 30°-40°to line-of- 
sight). From equation (|4)), we find Doppler factor is a slowly 
varying function of the magnetic field, i.e., 6 ~ _B(1)^''^^. 
Then the uncomfortable small viewing angles of the outer 
jet may be the artifact of the assumed equipartition state 
of the magnetic field down the whole outfiow. M oreover, we 
notice the research o fH emz fc BegehnanI |[l993) on M87 jet 
suggests that there is some departure from equipartition for 
the magnetic field of knot A, i.e., 0.2 < B/Bcq < 0.6. Thus 
we can further evaluate the effect of magnetic field values 
on the derived jet orientation based upon their constraint 
to knot A. 

W ith the constraint to knot A (|Heinz fc BegelmanI 
Il997h . we assume the magnetic field of this jet region with 
no beaming, -B(l), has some departure from the related 
equipartition value (see column 5 of Table [2} . Thus we have 

B(l)=CBcq(l),(0.2<e<0.6) (7) 

where ^ is a dimensionless parameter. Using different val- 
ues of 5, or different -B(l), we then repeat the same proce- 
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Figure 3. Plots of (a) observed apparent speed (in unit c), (b) Doppler beaming factor, (c) intrinsic bulk Lorentz factor and (d) viewing 
angle to line-of-sight versus projected distance down M87 jet. 



dures described above to calculate the related Doppler fac- 
tors, bulk Lorentz factors and viewing angles. The results 
are listed in Table |4l 

Inspection of Table |4] suggests that the orientation of 
knot A depends on the level of departure from equipar- 
tition. If 0.4 < ^ < 0.6 (small departure), its orienta- 
tion can maintain the same viewing angle with the in- 
ner jet, i.e. 9 ~ 10° — 20°. Meanwhile, larger departure 
from cquipartition (0.25 < ^ < 0.3) leads to large incli- 
nation (9 ^ 30° — 40°) co nsistent with that inferred by 
feicknell fc BegelmanI l|l996l ). Based on these results, we ar- 
gue that the outer jet of M87 starting from knot A suffers 
some departure from equipartition. This departure may re- 
sult from accumulative mass entrainment during jet propa- 
gation. 



5 CONCLUSIONS AND SUMMARY 

Based on multi-waveband fitting with moderate synchrotron 
spectrum models and radiation-escape balance at the break 
frequency, we derive the distribution of Doppler beaming 
factors down M87 jet. Coupling the data of proper motions, 
the intrinsic velocity and orientation are calculated for eight 
knotty regions, which presents some hints for the overall 
dynamics of this famous extragalactic jet. The main results 
are summarized as following: 

(1) The derived distribution of Doppler beaming factors 
gives good support to the 'modest beaming' scenario, with 
5 varying between ~ 2 — 5. 

(2) The bulk flow of M87 jet reveals sharp deceleration 



at the inner jet, which may be caused by continuous en- 
trainment of external gaseous medium during propagation. 
Meanwhile, the intrinsic speed of the outer jet remains 
roughly constant after early deceleration. 

(3) The viewing angles of the inner jet regions are fully 
consistent with 10°-19°to line-of-sight suggested by previous 
HST proper motions study of M87 jet. The outer jet, how- 
ever, shows systematic deflection off the inner jet to much 
smaller inclination (9 <^ 10°). Further calculation of knot A 
suggests this deflection can be regarded as evidence that the 
outer jet suffers some departure from equipartition. 

(4) The nucleus region of M87 jet should have an viewing 
angle close to its first knot HST-l, i.e. 9 ~ 15°. Then our 
result favors the idea that M87 may be a misaligned blazar. 



We 



should 



also 



comment 
Using a 



the recent work of 
relativistic jet model. 



Laing fc B ridle (200. 

Laing fc Brid le ( 20081 ) try to model the high resolution VLA 
images of five FR I radio galaxies other than M87. They find 
that a decelerating jet model can fit the observed bright- 
ness and polarization structures well. Moreover, they suggest 
the deceleration is most likely caused by mass entrainment. 
Their results provide strong support to our result, i.e. the 
decelerating velocity profile found in typical FR I jet of M87 
is reliable. 

Once the related jet parameters are derived in Table ([5} 
and Table ([3]), we can estimate the contribution to jet power 
carried by the radiat ing electrons and the magnetic field as 
(|Kataoka et al.ll2008l ) 



Lj«2^i?^cr^(|-) 



(8) 
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Table 4. Kinetic Parameters of Knot A with Different Magnetic Field 





5 


r 


e (deg) 


0.25 


1.42 ±0.57 


1.15 ±0.11 


38.6 ±35.5 


0.3 


1.58 ±0.63 


1.19 ±0.16 


30.1 ± 28.2 


0.4 


1.87 ±0.75 


1.27 ±0.24 


20.3 ± 19.0 


0.5 


2.13 ± 0.85 


1.36 ±0.31 


15.0 ± 13.7 


0.6 


2.36 ±0.95 


1.45 ±0.37 


11.9 ± 10.7 



where the factor 2 means the total internal energy density is 
twice of the energy density of equipartition magnetic field. 
Making use of derived parameters of HST-1 knot and the 
transformation formulae in Appendix A, we get Lj ~ 9 x 
10''^erg.s~^. Our result agrees well with Lj « lO^'^erg.s"^ 
suggested by different methods of previous rese arches (e.g., 
iBicknell fc BegelmanlFigQel . IStawarz et al.ll2006l '). 

Another issue we should confirm is whether our derived 
jet parameters are consistent with the conservation of the 
magnetic flux along M87 jet. In the jet frame, the mag- 
netic flux 5* = B'{nR'^). Using derived jet parameters in 
Table ([2]) and transform with formulae of Appendix A, we 
can estimate the related magnetic fluxes (in /iG.pc^) as (the 
equipartition level of B in outer jet is somewhat uncertain, 
thus we only give results for inner jet regions): 7.23 x 10^ 
(Nucleus), 2.14x10^ (HST-1), 4.10 x 10^ (D-east), 5.63 x 10^ 
(E), and 1.07 x 10® (F). The results indicate the magnetic 
fluxes remain roughly conservative down the jet. 

Finally, we should test the validity of main model as- 
sumption on the dominance of synchrotron cooling process. 
To compare the dominance of synchrotron cooling rate over 
SSC (synchrotron self-Compton) cooling rate, we have the 
ratio 

(d7/dt)ggc ^ 

where Ug — B'^/Stt is the intrinsic energy density of mag- 
netic fleld. The local energy density of photon field, f/pj^, can 
be estimated with intrinsic total synchrotron luminosity of 
the jet region, i.e. U^i^ = L'g/{4:nR'^c). Making application 
to the jet parameters and observed total luminosity of jet 
regions, and transform with Appendix A, we find out the 
ratio defined in equation @ range between 10~^ — 10"'^ 
along M87 jet. The results are consistent with the model 
assumption on the dominance of synchrotron cooling. 

Meanwhile, for the inner parts of the jet, the M87 host 
galaxy may provide relative intense photon fleld for the IC 
scattering process, resulting in an additional radiative cool- 
ing of the ultra-relativistic electrons. We should evaluate 
this possible effect in details. The star light from M87 host 
galaxy contributes to the jet as an incident external soft 
photon field at the optical band. To simplify the analysis, 
we assume the external photons have some characteristic 
frequency fext « 5 x 10^* Ifz. Then after external IC scat- 
tering (EC process) by ultra-relativistic electrons in the jet, 
the emerged high- energy photon s will have some character- 
istic frequency as (jFossatill 19981 ) 

l/EC ==; ^T^'^Fl^ext, (10) 



where 7 is the Lorentz factor of the ultra-relativistic elec- 
tron, S and r refer to Doppler beaming factor and bulk 
Lorentz factor of the jet flow, respectively. 

As to the ultra-relativistic electrons responsible for syn- 
chrotron X-ray emissions in large-scale jets, their typical en- 
ergy should be 7 ~ 10® or even higher. According to our 
inferred jet parameters of M87, we find 5T ~ 10 for the in- 
ner jet regions. Therefore, the resulting high-energy photons 
from EC process locate around i/ec ^ 10^^ Hz, which falls in 
the TeV energy band. The result above suggests that signif- 
icant EC cooling of ultra-relativistic electrons from the host 
galaxy will predict significant contribution to t he observed 
TeV s pectra of M87. However, detailed study of lLiu fc Sheiil 
(|2003) implies that the large-scale knots are unlikely to be 
the site for the TeV emissions recently dete cted in M87. 
Through modelling of the core jet of M87, E enain et al.l 
(2008) do not also find evidence of signiflcant contribution to 
TeV emissions from EC process of M87 host galaxy. There- 
fore, we argue that the synchrotron cooling dominated as- 
sumption of the model should hold in the inner jet regions 
of M87. 

During the preparation of th is paper, we notice the lat- 
est work of ISahavanathanI (|2008l ) on M87 jet knots. An im- 
proved two-zone (emission region and near by acceleration 
region) CI synchrotron model is pr esented b y Sahava nathanI 
(2008), very similar to the work of lLiu fc Shea (2003). Both 
of the modified CI models can fit the broad-band SED much 
better than standard synchrotron models for most of the 
M87 jet knots. We comment, however, another possibility 
for further improvement of the model. All of the current 
synchrotron models only consider synchrotron radiation of 
relativistic electrons spiralling in an uniform magnetic field, 
which is unlikely to be the real case. 

It is generally believed that large-scale jet knots 
are places of strong shocks, where ultrarelativistic elec- 
trons are accelerated via the first-order Fermi pro- 
cess (e.g., iBlandford fc Eichleilll987l : lOstrowski fc Bednar3 
[2002"). There is no strong evidence to support the sugges- 
tion that the magnetic field lines behind the shock front are 
straight and uniform. The structure of the magnetic fleld 
may be highly tangled, composed of curved magnetic fleld 
lines with different scales of radius. 

On the other hand, theoretical researches on 
synchro tron radiation in c urved magnetic fleld 
lines (IZhang fc Chengj Il995l . ICheng fc ZhaiiS Il996l . 
IZhang fc YuanI 1 19981 : the mechanism is named 'synchro- 
curvature radiation' by the authors) suggest that the 
curvature of magnetic fleld lines plays an important role 
in the spectral slope of the synchrotron spectrum, i.e. 
at the high-energy portion of the SED. Recently, such 
a new mechanism was used to explain successfully the 
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high-energy excess of several gamma-ray bursts (GRBs; 
iDeng. Xia fc Liul |2005| ) the broad-band emission of which 
may also originate from relativistic jets. Since standard 
synchrotron models show systematic deviation from the 
observed X-ray spectra of M87 jet knots, it would be inter- 
esting also to explore the effect of curvature of magnetic 
field lines on the broad-band SED of the M87 jet and other 
objects in possible future work. 
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APPENDIX A: RELATIVISTIC 
TRANSFORMATIONS TO DERIVE 
DOPPLER BEAMING FACTOR 

In this section, we perform relativistic transformations from 
the jet rest frame (primed) to the observer frame (unprimed) 
to derive Doppler beaming factor of the jet knot. Firstly, let 
us assume that the jet knot, observed as a spherical emit- 
ting region with the radius 7? and the deprojected observed 
length L ~ _R, is a moving blob. 

In the jet rest frame, the synchrotr on cooling time a t 
the break frequency can be expressed as l|Pacholczvklll970l ) 



isyn(^B) = 5.08 X 10^^ B' 



-3/2 /-1/2 



(Al) 



where the equipartition magnetic field B'^q, is in fiG, and 
the break frequency i/g, is in GHz. 

While performing Lorentz transformation from jet 
frame (primed) to the observer frame (unprimed), the 
break frequency transforms as (see formula (C6) of 
iBegelman. Blandford fc Reej|l984 ) 

J^B = I^B,ohs/S. (A2) 

Meanwhile, we use equation (A8) of IStawarz et al.) (|2003l ) 
to transform the equipartition value of the magnetic field, 
which corrects the formula (A7) of [Harris fc Krawczvnskil 
l|2002l ). The transformation is expressed as 



Beq = B{1)5 



-5/7 



(A3) 



where B{1) is the equipartition magnetic field calculated for 
no beaming {5 — 1). 

We then substitute equations (|A2|I and HA3|) into (|A1|) . 
and the synchrotron cooling time is rewritten as 



= 5.1 



X W^'^Bil) 



-3/2 -1/2 cll/7 



(A4) 

On the other hand, in the jet rest frame, L' — L/S ~ 
R/S. Therefore we have the volume of the blob V' = 
ttR^L' = TvR^L/S = V/S. For the assumed geometry of the 
knot, the spatial scale of interest in the escape ti mescale is 
L' . T hen we rewrite the form of escape timescale (|Kataokal 
1 19991 ) in the jet rest frame as 



tesc = vL' /c, 



(A5) 



(A6) 



where 77 is a dimensionless parameter, and c is speed of light 
in cm • s~^. 

The relation above then gives 

, L R 

As to M87 jet, we have FWHM measurement of the 
major axis (6^3) and minor axis (Sb) for the jet knots 
from high resolution im aging observations by Chandra 
(IPerlman fc Wilson|[2005h . Therefore, we can use the quan- 
tity y/9a, X 9h (in arcsec) to estimate the mean angular di- 
ameter of the region, which in turn gives average observed 
radius -Robs with the known distance of the object. The pro- 
jection scale of M87 jet is 1"— 77.6 pc, thus 



Ro 



0.5\/6'a X X 77.6 X 3.09 x 10^** cm. 



(A7) 



At the break frequency of the synchrotron spectrum, the 
synchrotron cooling time equals to the e scape time of elec- 
trons from the main emis sion blob (e.g.. llnoue fc Takaharal 
1 19961 . iKataoka et aklbOOOl ) . When we equal the right side of 



equation (|A4|) and (|A6[) . and rewrite vs.ohs in Hz, the final 
formula of S can be easily expressed as 



1.61 X 1021 



(A8) 



This paper has been typeset from a Tp^X/ ET]eX file prepared 
by the author. 



